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The adhesion bonds connecting a lipid bilayer to an underlying surface may undergo a condensa¬ 
tion transition resulting from an interplay between a short range attractive potential between them, 
and a long range fluctuation-induced potential of mean force. Here, we use computer simulations 
of a coarse-grained molecular model of supported lipid bilayers to study this transition in confined 
membranes, and in membranes subjected to a non-vanishing surface tension. Our results show that 
confinement may alter significantly the condensation transition of the adhesion bonds, whereas the 
application of surface tension has a very minor effect on it. We also investigate domain formation in 
membranes under negative tension which, in free membranes, causes enhancement of the amplitude 
of the membrane thermal undulations. Our results indicate that in supported membranes, this 
effect of a negative surface tension on the fluctuation spectrum is largely eliminated by the pressure 
resulting from the mixing entropy of the adhesion bonds. 

PACS numbers: 


I. INTRODUCTION 

Lipid bilayers serve as a physical barrier separating the 
content of the cell from its surroundings. They are em¬ 
bedded with different types of proteins responsible for 
numerous biological functions. An important class of 
membrane proteins are cell adhesion molecules that form 
specific (ligand-receptor) bonds with different biological 
elements, e.g., the extracellular matrix, the cytoskeleton, 
and neighboring cell membranes [l| . Cell adhesion plays 
a central role in many biological processes, such as T- 
cell activation as part of the immune response [H, cell 
migration Q and tissue formation Q . 

When membranes adhere to another surface, the ad¬ 
hesion bonds may condense to form large adhesion do¬ 
mains which not only provide mechanical stabil¬ 

ity but also promote cellular signalingpathways neces¬ 
sary for many biological processes [8|,l9|. The aggrega¬ 
tion of adhesion bonds is driven by attractive forces, in¬ 
cluding electrostatic and Van der Walls forces [l 0 [, inter¬ 
actions resulting from cytoskeleton remodeling [^, and 
membrane-mediated interactions M- The latter mech¬ 
anism stems from the suppression of membrane ther¬ 
mal fluctuation caused by the adhesion. The associated 
decrease in fluctuation entropy is minimized when the 
bonds cluster together into a single adhesion domain. 

Over the past two decades, considerable efforts have 
been made to better understand the role of membrane 
mediated interactions in the formation of adhesion do¬ 
mains. To this end, a common theoretical approach is the 
use of discrete lattice models, where the occupied sites 
represent membrane adhesion bonds, while empty sites 
stand for free fluctuating membrane segments |13l-[l^. 
Computer simulations of such models have shown that 
membrane fluctuations can facilitate conditions for con¬ 
densation of adhesion bonds. Nevertheless, it has been 
also concluded that the fluctuation induced interactions 


alone are too weak to promote aggregation and, hence, 
other attractive interactions between the bonds must also 
exist. A similar conclusion was recently reached in sim¬ 
ulation studies employing coarse-grained (CG) molecu¬ 
lar models, which provide a more physical description of 
the system. In one such study, we simulated a bilayer 
membrane with a small fraction of lipids from the lower 
monolayer connected to an underlying surface, and addi¬ 
tionally introduced a short range attractive potential of 
depth e between the adhesive lipids [ 2 ^ . Our simulations 
revealed that the adhesive lipids underwent a first order 
condensation transition when the strength of the short 
range potential between them exceeds a threshold value 
ec > 0. The fact that adhesion domains do not form 
for e = 0 implies that the fluctuation entropy gained by 
the aggregation of the adhesive lipids only partially com¬ 
pensates for their loss of mixing entropy. In another GG 
simulation study, Noguchi demonstrated that cell junc¬ 
tions connecting more than two membranes can aggre¬ 
gate without an additional attractive potential between 
them HH. This observation can be understood because 
the total entropy of thermal undulations increases lin¬ 
early with the number of membranes, while the mix¬ 
ing entropy of the cell junctions is independent of 

Previous simulation studies have focused on the aggre¬ 
gation behavior of adhesion bonds in tensionless mem¬ 
branes. Here, we wish extend our investigations to lipid 
bilayers subjected to surface tension and confinement. 
Under such conditions, the membrane long wavelength 
thermal undulations are suppressed, which weakens the 
fluctuation induced attraction. Specifically, a harmonic 
confining potential of strength 7 > 0 suppresses thermal 
fluctuations at length scales larger than ~ (k/ 7)^/^, 
where k is the bending modulus of the membrane. Like¬ 
wise, surface tension r > 0 affects bending modes with 
wavelengths greater than ^ The influ¬ 

ence of tension and confinement on the fluctuation me¬ 
diated potential of mean force (PMF) was analyzed for 
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the case of two adhesion bonds [ 22 |. It was shown that 
the PMF depends logarithmically on the pair distance: 
0 (r) = CkBThi{r)^ where C = 2 at short distances, 
C = 1 at long distances [r ^ ^r) for stressed membranes, 
and C = 0 for long distances (r ^ for confined 
membranes. These findings suggest that tension and 
confinement may hinder the formation of adhesion do¬ 
mains, and that the phase transition into the condensed 
state requires a stronger short range attraction between 
the bonds. Therefore, one should expect the transition 
threshold value Cc to grow with r and 7 ; however, the 
magnitude of the increase in ec remains hard to assess 
due to the many body nature of the PMF 0 - One may 
also consider the case of a bilayer membrane under neg¬ 
ative tension. Negative tensions amplify the amplitude 
of the long wavelength undulations (and, for very large 
negative values, even lead to instability), and this should 
cause a reduction in ec- It is especially interesting to ex¬ 
amine whether Cc can decrease to zero, in which case the 
adhesion domains will form without an additional short 
range potential, i.e., on purely entropic grounds. As will 
be shown further on, the results of the current study sug¬ 
gest that confinement might have a very strong impact 
on ec, while the application of a positive surface tension 
hardly influences the transition into an aggregated phase. 
In cases where the membrane is subjected to a negative 
surface tension, we observe the formation of elongated 
adhesion domains close to the transition point. 


II. MONTE CARLO SIMULATIONS 


In order to address the above issues, we conducted 
Monte Carlo simulations of a lipid bilayer using the 
Cooke-Deserno implicit-solvent CG model 1^. The de¬ 
tails of the simulations can be found in ref. [^. Briefly, 
we simulated bilayers of 2N = 2000 lipids (1000 lipids per 
monolayer), where each lipid consists of one hydrophilic 
(head) bead and two hydrophobic (tail) beads of size a. 
A flat impenetrable surface was placed below the mem¬ 
brane at z = 0 , and adhesion was established by con¬ 
necting Nf) head beads from the lower monolayer to it. A 
Lennard-Jones type potential with a tunable depth e was 
introduced between the adhesive head beads. We simu¬ 
lated systems with various densities of adhesive lipids, 
(j) = Ni)/N^ and for different values of the interaction po¬ 
tential depth e. Our simulations for confined membranes 
included an additional impermeable surface placed above 
the upper monolayer at ^ = ^^conf- Simulations of mem¬ 
branes under constant mechanical surface tension r were 
carried out according to the method described in ref. [ij . 

Typical equilibrium configurations of tensionless mem¬ 
branes with (j) = 0.2 are shown in Fig. [T]for e = 0.2/cbT 
(A) and e = lAksT (B). In the former, the adhesive 
lipids (whose head beads are shown in the figure in red) 
are scattered across the surface, while in the latter they 
condense into a large cluster. Fig. [IJC) shows the av¬ 
erage energy E of pair interactions between adhesive 



FIG. 1: Bottom view of equilibrium snapshots of a tensionless 
bilayer with 0 = 0.2 showing (A) a scattered distribution of 
adhesion bonds for e = OAksT, and (B) a condensed phase for 
e = lAksT. The head and tail beads are colored in white and 
blue, respectively, while the adhesive beads are colored in red. 
(C) The average number of contacts between the adhesion 
bonds as function of e. The condensation transition value Cc 
is empirically defined via the equality (Nc) — 1.5. The line 
serves as a guide to the reader’s eye. 


beads, normalized per bond and expressed in dimen¬ 
sionless units by dividing by the potential strength e. 
This quantity, which is a measure for the typical number 
of contacts between the adhesive beads (Ac), increases 
sharply when the condensation transition takes place. 
From the results of Fig. [HC), we empirically determine 
the transition value as the one for which {Nc) ^1.5. 


III. RESULTS 

A. Membranes under confinement 

We first start with our simulation results for non 
stressed membranes confined between the supporting un¬ 
derlying surface (on which the adhesive beads reside) 
and a second upper surface. The former is located at 
z = 0 , just underneath the tips of the head beads of the 
lower leaflet, while the latter is placed at z = ^^conf ^ ba. 
The degree of confinement increases when Zconi decreases 
which, in turn, would lead to stronger suppression of 
the membrane thermal fluctuations and a shift in the 
transition threshold Cc to larger values. This trend is 
demonstrated in Fig.[2l showing our simulation results for 
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FIG. 2: The average number of contacts between the adhesion 
bonds in membranes confined by a surface located at ^conf = 
9cr (green diamonds), ^conf = 7.5a (blue triangles) and ^conf = 
6a (red squares), as a function of e, for (A) 0 = 0.1, and (B) 
0 = 0.2. Results for non-confined membranes are denoted 
by black circles. The lines serve as a guide to the reader’s 
eye. The statistical errors are comparable to the size of the 
symbols. 


Zconi = 6(7, 7.5(7, and 9(7 in supported membranes with 
0 = 0.1 (A) and (j) = 0.2 (B). For ^^conf = 9(7, our results 
(green diamonds) match perfectly with the results ob¬ 
tained for non-confined membranes (black circles). This 
means that the rate of collisions between the membrane 
and the upper surface is negligibly small and, therefore, 
it does not affect the fluctuation spectrum. Lowering 
the confining surface by a distance equal to the size of a 
bead and a half to Zconi = 7.5(7 has a more noticeable ef¬ 
fect on membrane thermal undulations, which leads to 
an increase in the condensation transition value from 
Cc — O.GbkBT for non-confined membranes at (j) = 0.1 
to Cc — O.SksT. For 0 = 0 . 2 , the shift is smaller, from 
Cc — O.GksT to Cc — O.TksT. When the upper surface is 
further lowered to ^^conf = 6(7, it touches the tips of the 
head beads in the upper leaflet, as the thickness of the 
bilayer is equal to the size of six beads. A confining sur¬ 
face located at Zconf = 6(7 completely suppresses thermal 
undulations, and eliminates the fluctuation mediated in¬ 
teractions between the adhesion bonds. Under these con¬ 
ditions, the threshold value for aggregation increases to 
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FIG. 3: The average number of contacts between the adhesion 
bonds in membranes under mechanical surface tension r = 
0.48 A:sT/( 7^, as a function of e for (A) 0 = 0.1, and (B) 
0 = 0.2. The results are plotted in red squares, and are 
compared with the results for tensionless membranes (r = 0) 
that are depicted by black open circles. The lines serve as 
guides to the eye. The statistical errors are comparable to 
the size of the symbols. 


Cc — 1 . 2 kBT at 0 = 0.1, and Cc — IkBT for (j) = 0.2. 
These values are approximately twice larger than the 
corresponding values found when no upper plate exists 
(cc — O.GbkBT and Cc — 0 . 6 kBT, respectively), which is 
in accord with the conclusions of refs. [I^, [^, that the 
entropic gain of aggregation compensates for, roughly, 
half of the loss in mixing entropy of the adhesion bonds. 


B. Membranes under surface tension 

We now proceed to our simulation results for non- 
confined membranes subjected to a lateral surface ten¬ 
sion r > 0. Since, just like steric confinement, surface 
tension also suppresses the thermal undulations of the 
membrane, one may expect to also find an increase in the 
value of €c. However, our computational results that are 
depicted in Fig[ 3 l demonstrate that the (Nc) vs. e curve 
remains almost unchanged when the tension increases 
from r = 0 (black open circles) to r = 0 . 48 /cbF/( 7 ^ (red 
squares) for both 0 = 0.1 (A) and (j) = 0.2 (B). This ob- 
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FIG. 4: The average number of contacts between the adhesion 
bonds in membranes under surface tension r = — 0.24A:BT/cr^, 
as a function of e for (A) 0 = 0.1, and (B) 0 = 0.2. The 
results are plotted in red squares, and are compared with the 
results for tensionless membranes (r = 0) that are depicted 
by black open circles. (C) The mean projected area per lipid 
as a function of e, for 0 = 0.1 (black) and 0 = 0.2 (red). 
Circles denote the results for tensionless membranes, while 
the results for r = — 0.24A:sT/cr^ are shown in squares. The 
lines serve as guides to the eye. 


servation suggests that the fluctuation spectrum of the 
simulated membranes is very little affected by the ten¬ 
sion, which is expected given that the bending modulus 
of the model membrane has been previously estimated 
to be /i: SksT [25[. For r = 0ASkBT/a‘^^ the wave¬ 

length of the modes that are suppressed are larger than 


— 25(7, which is much larger than the 
typical distance d = 4(7 between the ad¬ 

hesion bonds when the latter are uniformly distributed 
on the surface (where Ap is the projected area of the 
membrane). In order for the effect of tension to be ap¬ 
preciable, one must apply stronger tensions; however, the 
application of a stronger tension would open pores in the 
membrane and tear it apart. In simulations with a sur¬ 
face tension of r = 0.64/c^T/(7^, we observed pores in 
most of the simulated membranes, and in the few that 
had managed to survive the tension, we still did not ob¬ 
serve any considerable changes in the strength of mem¬ 
brane mediated interactions. 

We next aimed to address the implications of applying 
a negative surface tension. In general, a negative tension 
imposed on the membrane leads to a reduction in its pro¬ 
jected area, Ap, and amplifies the long wavelength bend¬ 
ing modes. Hence, the fluctuation-induced attraction 
between the adhesion bonds is expected to be stronger 
which, in turn, implies that the threshold for condensa¬ 
tion Cc should become smaller than in tensionless mem¬ 
branes. To test this hypothesis, we simulated the mem¬ 
brane under a negative tension of r = —0.24/cpT/(7^. We 
performed two sets of independent MC simulations, one 
starting from a random distribution of adhesion bonds, 
and another where initially the adhesion bonds were or¬ 
ganized in one large cluster. The system was equilli- 
brated until configurations originating from these two 
distinct initial conditions achieved similar characteris¬ 
tics. Fig. |4] shows our results for {Nc) as a function 
of e for 0 = 0.1 (A) and (j) = 0.2 (B). Contrary to our 
expectation to observe a reduction in ec, the data for 
r = —0.24/cpT/(7^ appears almost identical to the results 
of the tensionless case, with Cc — 0.6/cpT for both values 
of (j). Taken together with the data for r = -h0.48/cpT/(7^ 
(see Fig. [3|), we conclude that application of surface ten¬ 
sion has a negligible effect on the condensation transition, 
within the range of tensions where the supported mem¬ 
brane is mechanically stable. 

The negative tension, however, does have an impact on 
the shape of membranes. Freely fluctuating bilayers as¬ 
sume buckled configurations at negative tensions larger 
(in absolute value) than Tc — —^-'K^njAp [26|, In 

this study, ApjN 1.33(7^ [see Fig. |11[C)], which gives 
Tc — —0.24/cpT/(7^. In ref. [^, a similar model mem¬ 
brane consisting of the same number of lipids was sim¬ 
ulated and, indeed, for r = Tc the membrane appeared 
quite buckled. In supported membranes, however, the 
emergence of buckled configurations occurs only in mem¬ 
branes with large adhesion domains (i.e., for e > ec). 
Fig.[5]shows typical equilibrium configuration for (j) = 0.1 
with e = 0 (A), 0.6/cpT (B), and 1.0/cpT (C). Each con¬ 
figuration is shown both in side and bottom views. When 
the attractive potential is set to e = 0, the distribu¬ 
tion of the adhesion bonds is scattered and the mem¬ 
brane remains fairly fiat. This indicates that the mix¬ 
ing entropy of the bonds dominates the fluctuation en¬ 
tropy of the bilayer, despite the imposed negative ten- 
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(A) (B) (C) 



e = 0 e = O.G/cbT e = l.O/c^T 


FIG. 5: Typical equilibrium configurations of membranes un¬ 
der surface tension r = —0.24 A;bT/( 7^ with density of adhe¬ 
sion bonds <f> = 0.1 for (A) e = 0, (B) e = O.G/csT and (C) 
e = l.^ksT. The figures in the upper and lower rows dis¬ 
play bottom and side views of the system, respectively. Color 
coding is the same as in Fig. [H 


sion [see Fig. EfA)]. For e = the short range 

pair interactions between the bonds lead to their aggre¬ 
gation. Once the bonds condense, their influence on the 
thermal behavior of the membrane is greatly weakened, 
and strong bending undulations appear [see Fig. [5][C)]. 
Close to the condensation transition, at e = 0.6 /cbT, the 
system exhibits some interesting features: The ampli¬ 
tude of one of the two longest wavelength bending modes 
[with wavevector ^(i,o) = (27r/\/Ap)(l, 0), or ^(o,i) = 
(27r/A/Ap)(0,1)] grows considerably, and the membrane 
assumes an anisotropic buckled configuration. The ad¬ 
hesion bonds are concentrated throughout the minimum 
of the dominating bending mode, forming an elongated 
domain (“stripe”) [see Fig.[5][B)]. These observed charac¬ 
teristics represent an intricate balance between the driv¬ 
ing forces that govern the thermodynamic behavior of the 
system. Under negative tension, the system benefits from 
a reduction in the projected area, leading to a decrease in 
the Gibbs free energy. The membrane, however, is quite 
incompressible and, thus, the reduction in Ap must be 
accompanied by an increase in the area stored in thermal 
fluctuations whose amplitudes grow. The modes that ex¬ 
perience the largest increase in amplitude are the softest 
ones, corresponding to g(i,o) and q(o,i) [HI (see footnote 
ref. [29|). In our simulations, we seldom observed situ¬ 
ations where both these modes were simultaneously ex¬ 
cited [ 3 ^, which can be linked to the mixing entropy of 
the adhesion bonds. When only one of the long modes is 
dominant, the contact area between the membrane and 
the surface, which is available for the presence of the ad¬ 
hesion bonds, is larger than in configurations where both 
modes are excited. 

The interplay between the mixing entropy of the adhe¬ 
sion bonds and the contribution of the negative tension to 
the free energy is further demonstrated in Fig. |4l[C), de¬ 
picting the mean projected area {Ap) (normalized by the 


(A) (B) 



FIG. 6: Gonfigurations of membranes with 0 = 0.1 and e = 
1.0/cpT under a strong negative tension r = —0.32/cpT/cr^, 
showing a spherical protrusion (A) and a tubular one (B). The 
upper and lower rows display top and side views of the mem¬ 
brane respectively. The tail and adhesive beads are colored 
in blue and red, respectively. The head beads are colored in 
grayscale to reflect their height above the surface, with lighter 
colors representing a higher bead. 


number of lipids per monolayer N = 1000) as a function 
of e for r = —0.24/cpT/<j^ (squares) and r = 0 (circles). 
In the tensionless case, we observe a very mild decrease 
in (Ap) as e increases, occurring mainly around Cc- For 
r = — 0.24/cpT/cr^, (Ap) maintains a value close to the 
tensionless case for e < Cc, and drops significantly for 
e > Cc. In the latter regime, we also observe an increase 
in the area fluctuations, resulting in larger uncertainties 
(error bars) in our estimates of (A^). The sharp decrease 
in (Ap) and the concurrent increase in the area fluctua¬ 
tions are anticipated outcomes of a negative surface ten¬ 
sion [ 2 ^. The fact that they are observed only above 
the condensation transition is consistent with the picture 
discussed in the previous paragraph that, below ec, the 
effect of the negative tension is largely eliminated by the 
pressure resulting from the mixing entropy of the adhe¬ 
sion bonds. Notice that the sharp decrease in the mean 
projected area and the increase in the area fluctuations 
are much more noticeable for (j) = 0.1 than for (j) = 0.2. 
This is to be expected because the smaller 0, the smaller 
the restrictions imposed by the adhesion bonds on large 
thermal undulations which are directly coupled to the 
projected area by the highly incompressible character of 
the membrane. 

Applying an even stronger negative tension causes the 
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membrane to lose its mechanical stability. This is demon¬ 
strated in Fig.[6l showing snapshots in top and side views 
of membranes with 0 = 0.1 and e = l.^ksT subjected to 
r = —0.32kBT/a^. In these snapshots, the head beads 
are colored in grayscale, with lighter colors indicating 
beads located higher above the underlying surface. The 
application of a strong negative tension causes the sup¬ 
ported membrane to develop large protrusions with either 
spherical [Fig. [6][A)] of tubular [Fig. [6l[B)] shapes. The 
adhesion bonds (which are colored in red, and are only 
partially visible) are concentrated in the periphery of the 
protrusion, where the membrane is in contact with the 
underlying surface. We note that the observed protru¬ 
sions tend to evolve slowly and, therefore, the snapshots 
shown in Fig. [6] may not represent true equilibrium struc¬ 
tures. On the other hand, we also note that very similar 
equilibrium structures, featuring spherical and tubular 
protrusions, have been recently observed in an experi¬ 
mental study where supported lipid bilayers were sub¬ 
jected to compression |^ . 

IV. SUMMARY 

The thermal fluctuations of a supported membrane in¬ 
duce an attractive PMF between the adhesion bonds, 
that promotes the formation of adhesion domains. In 
previous studies of systems where the adhesion bonds in¬ 
teract with each other via a pairwise short range attrac¬ 
tive potential, it has been concluded that the fluctuation- 
induced PMF significantly decreases the condensation 
transition point, ec, of the adhesion bonds. In this work, 
we investigate adhesion domain formation in stressed and 


confined membranes by using MC simulations of a CG 
supported bilayer model. Our simulations reveal that 
placing a plate above the membrane (in addition to the 
surface supporting the membrane from below), signifi¬ 
cantly alters the condensation point. In the most ex¬ 
treme case, where the membrane is completely flattened 
between the two plates, the value of ec increases by al¬ 
most a factor of two. In contrast to the impact of steric 
confinement, the application of a positive surface ten¬ 
sion, which also acts to suppress membrane undulations, 
hardly affects the transition within the range of tensions 
where the supported membrane is mechanically stable. 
Similarly, we find that a negative tension also has a very 
minor effect on the condensation transition. Neverthe¬ 
less, once the adhesion bonds are aggregated into a large 
domain (i.e., for e > ec), the negative tension affects the 
shape of the membrane and causes it to buckle. Below Cc 
the adhesion bonds are scattered across the membrane, 
which prevents the formation of buckled configurations. 
From the fact that ec is hardly affected by the applica¬ 
tion of a negative tension, we conclude that the mixing 
entropy of the adhesion bonds dominates the demixing 
entropy of thermal fluctuations. Close to ec, we observed 
the formation of elongated adhesion stripes. Such con¬ 
figurations emerge from an interplay between the mixing 
entropy of the adhesion bonds, the short range resid¬ 
ual potential, and the applied negative tension. Finally, 
under a very strong negative tension, we observe tubu¬ 
lar and spherical structures protruding out of the mem¬ 
brane’s plane, which indicates that the system is at the 
onset of mechanical instability. 

This work was supported by the Israel Science Foun¬ 
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